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A massive foreground cluster lens changes the shapes (shear effect) and number 
density (magnification effect) of the faint background galaxy population. In this paper 
we investigate how the shear, magnification and combined information can be used to 
constrain cluster mass profiles in the weak lensing regime. We develop maximum likeli- 
Q\ ■ hood techniques to make quantitative predictions for each of the methods. Our analytic 
c-| ! results are checked against Monte Carlo simulations. 

In general, we find that the shear method is superior to the magnification method. 
However, the magnification information complements the shear information if the former 
has accurate external calibration. For the magnification method, we discuss the effects 
c3 ■ of random and systematic uncertainties in the background galaxy counts. 
> 

H ' 1 Introduction 

Weak gravitational lensing has been recognized as a powerful tool to investigate the 
mass and mass distribution of clusters of galaxies (e.g., Webster 1985; Tyson et al. 1990; 
Kochanek 1990). In their pioneering paper, Kaiser & Squires (1993) pointed out that 
the tidal gravitational field of the cluster, as measured from the distortion of image 
shapes of the faint background galaxy population, can be used to reconstruct the two- 
dimensional projected mass profiles of clusters, without referring to a parametrized mass 
model. This method, modified in various ways later (e.g., Kaiser 1995; Schneider 1995; 
Seitz & Schneider 1995; Bartelmann et al. 1996; Seitz & Schneider 1996; Squires & Kaiser 
1996; Seitz & Schneider 1997, 1998; Seitz, Schneider & Bartelmann 1998; Lombardi & 
Bertin 1998a, 1998b) has been applied to more than a dozen clusters up to now (e.g., 
Fahlman et al. 1994; Smail et al. 1995; Squires et al. 1996a, b; Seitz et al. 1996; Luppino 
& Kaiser 1997; Clowe et al. 1998; Hoekstra et al. 1998; Kaiser et al. 1999). One of the 
main results from these studies is that in many clusters the projected mass distribution 
closely follows the distribution of luminous cluster galaxies. For some of these clusters, 
very large mass-to-light ratios have been reported, e.g., for MS1224. 7+2007, this ratio 
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has been determined by two independent studies to be larger than ~ 8OO/1 in solar units 
(Fahlman et al. 1994; Fischer 1999). 

One of the main difficulties of these studies is the existence of the so-called mass- 
sheet degeneracy (Falco et al. 1985; Schneider & Seitz 1995) which states that the image 
shapes are unchanged if the surface mass distribution k{6) (defined in the usual way) is 
replaced by k{6) — > \k{6) + (1 — A). In particular, this transformation keeps the critical 
curves of the lens mapping fixed, and thus the location of the giant arcs and arclets. The 
choice of the constant A in actual mass reconstructions is largely arbitrary. Provided the 
observed data field is sufficiently large, one might argue that the surface mass density 
has decreased to a value close to zero far away from the cluster center, thereby fixing 
A. However, for those investigations where the available data field is quite small, the 
invariance transformation implies a substantial uncertainty in the mass estimates. 

The mass sheet degeneracy affects the magnification as n(0) — > n(6)/X 2 ; hence, 
if magnification information can be obtained, the degeneracy can be broken. Two ap- 
proaches to measuring the magnification have been presented in the literature: Broad- 
hurst et al. (1995) considered the effect of magnification on the local number counts of 
faint background galaxies, employing the so-called magnification bias effect. Since the 
number counts of very faint galaxies are flatter than the critical slope (defined below), 
magnification leads to a depletion of the local number counts near the cluster center. 
This effect has been observed in at least two clusters (Fort et al. 1997; Taylor et al. 
1998). Alternatively, Bartelmann & Narayan (1995) have suggested using the size of 
faint galaxies at fixed surface brightness as a measure for the magnification, making use 
of the fact that gravitational light deflection conserves the surface brightness. In addi- 
tion, further measures of the weak lensing magnification have been proposed, such as the 
change in the redshift distribution of galaxies in a flux-limited sample (Broadhurst et 
al. 1995), or the effect on the two-point angular correlation function of distant galaxies 
(Moessner et al. 1998). 

Whereas the assumption underlying the cluster mass reconstruction from image 
ellipticities (the shear method, in what follows), namely that the intrinsic orientation 
of source galaxies are randomly distributed - or that there is no preferred direction of 
alignment, is simple and predicts that the mean ellipticity of galaxy images in a small 
region of the cluster equals the reduced shear, the magnification techniques need an 
external calibration: for the first of the methods mentioned above, the number density 
of background galaxies in a 'blank field' as a function of limiting magnitude (and other 
selection criteria, such as surface brightness etc.) needs to be known, and for the second, 
the mean size of galaxies as a function of surface brightness is required. Uncertainties 
in this calibration information immediately translates into uncertainties in the derived 
magnification. In addition, naive estimates of the accuracy of shear and magnification 
measurements (see Sect. 2) show that the shear method appears to be superior. 

Nevertheless, Broadhurst (1999) claimed that the magnification information can, at 
least in some cases, yield a more accurate determination of the projected radial mass 
profiles of clusters. In this paper, we investigate the substance of this claim in more de- 
tail. In Sect. 2 we briefly summarize the shear and magnification methods, and provide a 
rough estimate of their respective accuracy in the determination of local lens parameters. 
Sect. 3 describes a maximum-likelihood method for the use of shear and magnification 
information to determine mass profiles. The ensemble-averages of the likelihood func- 
tions are derived in Sect. 4; later, these are used to determine the characteristic errors of 
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parameter estimates. In Sect. 5 we describe the numerical simulations used to substanti- 
ate our analytic results, and the lens models that are considered. We present the results 
from the likelihood analysis in Sect. 6; when the unlensed background number density is 
known, we find that the shapes of the likelihood contours are different in both methods. 
In particular, for a wide range of situations, the slope of the mass profile is ill-determined 
from the shear method, and better constrained by the magnification method. In Sect. 
6 we also consider the influence of uncertainty in the number density on the likelihood 
analysis. We summarize our results and conclude in Sect. 7. 

2 Shear and magnification methods 

We use standard lensing notation; i.e., k(0) denotes the dimensionless surface mass 
density of the deflector - which is well defined if we assume that all source galaxies 
are at the same distance, an assumption which is reasonable for low-redshift clusters 
(za 0.25), and for many purposes not a bad approximation even for medium-redshift 
lenses. The deflection potential ip(0) is related to k through a Poisson-like equation, 
V 2- i/> = 2k. The trace-less part of the Hessian of ip describes the tidal gravitational field, 
which we summarize in the complex shear 7 = 71 + 1^2 = (V , ,n — ^,22) /2 + iV',12, where 
indices separated by a comma denote partial derivatives with respect to the position 
on the sky. The magnification of an image is the inverse of the Jacobian determinant of 
the lens equation, fx(0) = [det „4(0)] _1 , with det^4 = (1 - k) 2 - \^\ 2 . 

The shear method is based on the transformation between source and image ellip- 
ticity due to the tidal field of the deflector. Let e s denote the intrinsic ellipticity of the 
source, and e the ellipticity of the observed image. Here, we use the ellipticity parameter 
as defined in, e.g., Seitz & Schneider (1997), such that e is a complex number whose 
modulus, in the case of elliptical isophotes with axis ratio r is |e| = (1 — r)/(l + r), and 
whose phase is twice the position angle of the major axis. The locally linearized lens 
equation yields the transformation between source and image ellipticity, 



where g = 7/(1 — «) is the complex reduced shear, and the asterisk denotes complex 
conjugation. The transformation (1) is valid only in the non-critical parts of the lens, 
i.e., where det^4 > 0, to which we will restrict the discussion here. Let p s e (e s ) d 2 e s be the 
probability that the source ellipticity lies within d 2 e s of e s , then the expectation value of 
the image ellipticity is 



as shown by Schramm & Kayser (1995) and Seitz & Schneider (1997). This result implies 
that the ellipticity of each galaxy image provides an unbiased estimate of the local reduced 
shear (it should be noted here that the reduced shear is unchanged under the mass-sheet 
invariance transformation). 

The magnification method based on the local number counts of background galaxies 
arises from the magnification bias (e.g., Canizares 1982) which states that the local 
cumulative number counts n(6;S) above a flux limit S are related to the unlensed counts 
n (S) by 




(2) 
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n{9-S) = -n (-) . (3) 



If we assume that the number counts follow (locally) a power law of the form n oc >S _/3 , 
then 

n(0)=n o ^~ 1 (4) 

at any fixed flux threshold. This implies that if the intrinsic counts are flatter than 
1, then the lensed counts will be reduced relative to the unlensed ones. This number 
depletion is the signature of lensing in the magnification method considered here. 

Let us briefly consider the signal-to-noise of a shear and a magnification detection by 
those respective methods. The noise in the former case is due to the intrinsic ellipticity 
distribution, whereas it is at least Poisson noise for the magnification method. If the 
source galaxies have a significant angular correlation, the resulting noise will be larger 
than the naive Poisson estimate used below; we consider this in Sect. 6.3. 

Thus, assume that in some (small) solid angle the shear is constant, and that there 
are iV 7 background galaxy images for which the ellipticity can be measured. The signal 
is fif, and the noise is a e / ^/iV 7 , where a e is the ellipticity dispersion. Thus, the signal-to- 
noise for a shear detection is 

s\ \g\ 



(5) 



N, 
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Similarly, assume that the magnification is constant over a solid angle and that, in the 
absence of lensing, galaxies would be expected in that region. Then, the magnification 
signal is \AN\ = — l|iV M , and the noise is y/N^, ignoring the change of the number 

density in the noise estimate. Thus, 

|) =|^- 1 -1| V / ^- (6) 

Specializing now to the case of weak lensing, i.e., k < 1, I7I <C 1, so that 7 ~ g, the 
magnification becomes \i xs l + 2«, and a first-order expansion yields \AN\ ss 2n\l — /3|iV M . 
Then, the ratio of the signal-to-noise estimates becomes 



(S/N) 7 | 7 | 1 



(S/N)„ k 2<r e \l-P\ 




(7) 



In a typical situation, the magnitudes of the shear and the surface mass density are 
similar (they are exactly the same for a singular isothermal sphere), and so the first 
factor in (7) is of order unity. The slope of the galaxy number counts is flattest for the 
faintest flux limits and for redder colors, and for these, (3 ^ 0.5. The variance of the 
intrinsic ellipticity is of order a e 0.2, and so the second factor in (7) is of order 5. The 
last factor can be expected to be smaller than unity, since the flux limit for counts can 
be fainter than that for which ellipticity measurements are still possible, but in order 
to make up for the second factor, the number density used for the counts would need 
to be 25 times higher than for the shear estimates, which at the assumed slope of the 
number counts would correspond to about 6 magnitudes! This of course is unrealistic - 
the ratio of over iV 7 is more typically of order 3 or 4 - and we therefore conclude 
from this simple analysis that the shear method is considerably more sensitive than the 
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magnification method based on number counts. The other magnification method based 
on images sizes yields a signal-to-noise ratio comparable to the shear method, but since 
this effect has not been observed up to now, we shall not consider this method in this 
paper. 

3 Likelihood functions for shear and magnification 

We shall now consider the following situation: in a solid angle (e.g., around a cluster), one 
observes galaxy images at positions 0j, 1 < % < iV M , down to a given flux threshold. 
In addition, one observes N 1 galaxy images at positions 1 < % < iV 7 , for which an 
ellipticity can be measured; the two sets of galaxy images can have objects in common. 
If we now assume that the mass distribution of the lens is described by a parametrized 
model with parameters 7^, 1 < i < M, what are the best-fitting parameters for the 
observation? 

We solve this problem by defining a likelihood function and then maximizing it. 
The likelihood is defined as the probability of obtaining the observables given the model 
parameters 7Tj, and is considered as the probability for the model parameters iii given 
the observables (see, e.g., Press et al. 1992 for a discussion on this point). 

Starting with the magnification effect, the probability of observing iV M galaxy images 
at positions 6i can be factorized in the probability of observing a total of galaxies, 
and in the probability that, given that there are iV M galaxies, they are located at 6i. For 
a lens model with parameters 7Tj, the expected number of galaxy images is 

(N,)=n,J d 2 9H0)f-\ (8) 

as follows from (4) ; we have used to denote the unlensed number density, to distinguish 
it from the number density of galaxies n 1 for which an ellipticity can be measured, and 
the integral extends over the solid angle of the observations. In practical applications, 
this area can be quite complicated since accurate photometry of very faint objects is 
impossible near brighter objects - like cluster galaxies - which have to be masked out 
(see Taylor et al. 1998 for an impressive illustration of this effect). The probability of 
finding galaxies in the data field follows a Poisson distribution (if, as we shall assume, 
the source galaxies are uncorrelated), and the probability that the i-th galaxy is located 
at &i is proportional to n(0j) oc \[i{6i)\^~ x . Then, 

with 

P(N; (N)) = ^ exp(- (N)) . (10) 

One can now see that the denominator in the product on eq. (9) is (iV M ) /n M , and that 
therefore the factor (N^)^ of the Poisson distribution is cancelled. In addition, the 
factors N^l and do not depend on the lens model and yield just an irrelevant mul- 
tiplicative constant in the likelihood; therefore we shall drop this factor. Then, the 
log-likelihood function becomes 
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V=-ln£ M = n M / d 2 0[//(0)]0-Ml-0)X>M0i)- (11) 

i=l 

The best fitting model is now obtained by minimizing with respect to the model 
parameters which are contained in the function n(0). 

Next we turn to the likelihood function for the shear estimate. For that, we need 
the probability that the observed ellipticity at position -di is e, given that the reduced 
shear is g {"&%). Assuming a Gaussian intrinsic ellipticity distribution, 

PA e) -^(l-exp(-lK)) (12) 



the observed probability distribution can be obtained from 



pM9)=Pl(e s (e\g)) 



d 



2 £ s 



d 2 t 



Pt(e s (e\g))- 



M 2 -i) 2 

|^*-1| 4 ' 



as in Geiger & Schneider (1998). The corresponding probability density p e (ei) for each 
lensed galaxy in the catalog is determined, and the likelihood and log-likelihood functions 
are given by 

£7 = n^(^lf W) S ^7 = - Yl ln Pe( e i\9(&i))- 

i=l i=l 

We shall follow a somewhat simpler road in the analytic treatment by assuming that 
the observed ellipticity distribution is a Gaussian; this is a reasonable assumption if the 
intrinsic ellipticity distribution is a Gaussian and if the reduced shear is not too close to 
unity. The expectation value of e is g, and the dispersion in the two directions parallel 
and perpendicular to the shear are the same (see Appendix of Geiger & Schneider 1999) 
and given by 



a 2 = 2a\ = 2o\ = 2rr(l - \g^ I d|e s | " 1 fV£ jf > ( 13 ) 



-(1-M 2 ) 2 r 

Jo 



|e s 


S Pli 


k S |) 


1 - 




e s 


2^2 



as will be shown in the Appendix. Provided \g\ is much smaller than unity, or that 
the intrinsic probability distribution does not extend significantly out to |e s | ~ 1, the 
denominator in the integral can be replaced by unity, for which the dispersion of the 
observed ellipticities becomes 

a^(l-\g\ 2 )a £ . (14) 

In the analytical treatment of Sect. 4, we shall use the approximation (13) for simplicity 
(for cr e = 0.2, the approximation for a differs from the true value by less than 5% for 
\g\ <1, and less than 1% for \g\ <0.4). The likelihood function for the shear method is 
then simply 



- _^ 1 ( \ei-9(*i)\ 2 \ 



and the corresponding log-likelihood function becomes, dropping irrelevant additive con- 
stants, 
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JV 7 

i=l 

Again, the best-fitting lens model is obtained by minimizing t 1 with respect to the model 
parameters, which enter g and thus also a. 

The combined shear and magnification likelihood is then obtained by just multiplying 
the respective likelihoods, or adding the log-likelihoods, 

4ot = ^ + V ( 16 ) 



-gWI 



+ 21n<r[y(^)] 



(15) 



4 Ensemble-averaged log-likelihood 

For a given set of observables, the likelihood functions defined in the previous section 
can be calculated; the minimum yields the best-fitting model parameters, and the width 
of the likelihood function yields the confidence region for the parameter estimate from 
this data set. In order to determine the characteristic errors of parameter estimates, we 
consider here the ensemble-averaged log-likelihood function. 

Given a lens model - hereafter called the true model, and characterized by subscripts 
't' - the ensemble-average of a quantity X is provided by 



(X) := P(N; (N)t) 



N=0 



X 



(17) 



where the first term averages over the probability to find N galaxies when the expec- 
tation value from the true model is (N) t , the second term is the integration over the 
probability that galaxy number i lies at 0^, and the final term integrates over the el- 
lipticity distribution. To see how this procedure works, we start with the log-likelihood 



function for the magnification method, i.e., setting X 



Since l„ does not depend on 



the ellipticity of the images, the final integration in (17) yields simply a factor one. The 
probability for the positions of the galaxies is 



PtiPi) 



/d30[/it(0)p-l 



Then, 



AT M =0 



II / ^°>ih(e,) 



f 



with 



(N,) t = n,J d 2 9 MO)?- 1 



(18) 



(19) 



(20) 



The first term in the log-likelihood function (11) depends neither on Oi nor on the ac- 
tual number of galaxies, and thus the averaging operator in (19) leaves this term 
unchanged. The second term of £ M consists of a sum over individual galaxy images; 
therefore, for the j-th galaxy image, only the term % = j in the spatial averaging in 
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(19) contributes, the other integrations just yield a factor one. Since there are iV M equal 
terms, one obtains 

oo NT ( 21 ) 

+ (i -/?) £ Wt) j d » % ( fl)1 M J d2 ° M*)]^ 1 • 

The sum over iV M can now be performed, and using (20) we obtain 

=n^J d 2 6 [Mf- 1 + n^l -0) J d 2 6 [^(O)f- 1 m/x(0) . (22) 
From the derivative of with respect to the parameter 7Tj, 

°m = (1 _ ^ -i. ?m { M9)] *-. _ , (23) 



one immediately sees that attains a minimum when the 7Tj attain their true values, 
i.e., when fx = fx t . 

If one assumes that the likelihood function behaves approximately like a Gaussian 
near its maximum, then the second partial derivatives of contain the information 
about the confidence region, and so we define 

where the values of the n derivatives are evaluated at the true parameter values. 

Next, we consider the ensemble average of the log-likelihood function for the shear. 
For that we have to specify the probability distribution for the image ellipticities for the 
true model. As mentioned before, we shall use a Gaussian with mean g^fl), the true 
reduced shear at position i9, and dispersion given by (13), again with g = g t , and denoted 
by at- Then, from (15) and (17), we see that, since l 1 is a sum over galaxy images, the 
integration operators work term by term, and so we get N 7 identical terms from the 
integrations. This then immediately allows one to perform the sum over iV 7 , as before, 
so that 

(3-1 



(£) = (n) tdH M*^ 1 



2 I ( \z-9^)?\ ( < 2 iM 



and (N 1 ) t is defined in analogy to (20), with n M replaced by n 7 . The e-integration is 
readily performed, and one obtains 

= n„ J d'« \» (*)]'- ( lgW - »jg' + "? W + 2M«)) • (26) 
By taking a first partial derivative of (£ 7 ) with respect to the parameter 7Tj, 
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d(£ 1 
dm. 



d 2 tf lvt(#)} 



(3-1 



dY 

dm. 



with 



dY 

dm 



dg 

dm. 



(g*-gt) 



2[|s-2t| 2 + K 2 -a 2 )] da 



dm 



(27) 



(28) 



we see that (£ 7 ) takes a minimum if n = n t , so that g = g t and therefore a = a t - 

From the ensemble-averaged log-likelihood function, we estimate the dispersion of 
the parameters one obtains from a single realization of the observables by using the fact 
that asymptotically (in the limit of infinitely many data points) the distribution of 2A£ 
tends to a Xm distribution with M being the number of model parameters. Since the 
number of galaxies is finite, the likelihood function is not Gaussian, and this estimate 
can only be approximate. We shall demonstrate with simulations that it provides a very 
good approximation indeed. 



5 Simulations and models 

The situation we shall consider is that of an axially-symmetric mass distribution for 
the cluster. We assume that the central part of the cluster is a strong lens where arcs 
and multiple images are found, and from their modelling the Einstein radius #e of the 
cluster is known. We then assume that the weak lensing analysis is performed in an 
annulus with inner and outer radii 6>i n and 6> ou t, with 6>i n > 6>e- If /3 < 1 (we shall assume 
throughout the paper that (3 = 0.5), the number density of galaxies in the annulus is 
always smaller than the unlensed number density no, since \x > 1 outside the critical curve 
(magnification theorem, see Schneider 1984). Hence, we can use the rejection method, 
as explained below. 

5.1 Description of the simulations 

Given a specific lens model, galaxy images are distributed in the following way: First, 
the expected number of galaxies in the annulus in the absence of lensing is calculated, 
(No) = no7r(0Q Ut — 0f n ), and a number Nq is drawn from a Poisson distribution with 
mean (No). Then, these Nq galaxies are provisionally distributed in angular radius, 
by drawing a random number £ uniformly distributed in [0, 1] , and assigning a radius of 
9 = v^t^out ~~ ^in) + ^?n- Then, a second uniform deviate r) G [0, 1] is drawn, and a galaxy 
at is put into the final 'catalog' only if |/t(#)]^ -1 > r\\ otherwise it is discarded. With 
this rejection method, the distribution of galaxies follows the distribution (18). For each 
galaxy in the 'catalog', a source ellipticity e s is drawn from a two-dimensional Gaussian 
probability density distribution given by equation (12), and the corresponding image 
ellipticity is obtained from (1). As a result, we obtain a catalog containing radial positions 
6i and image ellipticities e^. For such a catalog, and an appropriately parametrized family 
of lens models, the log-likelihood functions £ M , £ 7 and £tot can be minimized, resulting in 
the best-fitting parameters. 

The number density of galaxies for the shear and the magnification method can 
be different; throughout we assume that n 7 = 30arcmin -2 , and = 120arcmin -2 
(the former of these values is the typical number density used for shear analysis from 
deep ground-based data, the latter is similar to those used by Fort et al. 1997 and by 
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Broadhurst 1999). Furthermore, we shall assume that the Einstein radius of the cluster 
is #e = 0.5arcmin, as is approximately the case for the cluster CL0024+17 to which the 
magnification method has been applied (Fort et al. 1997; Broadhurst 1999). The resulting 
confidence regions of the parameter estimates depend solely on the ratios 6>i n /6>E, 6> ou t/#E, 
and on the products 9\n 1 and 9 2 E n^, and so can be easily scaled for different values of 
the Einstein radius. 



5.2 Lens models 

We shall consider two families of two-parameter lens models. The first (Family A) is 
characterized by a power-law profile outside the Einstein radius, i.e., 

k{9) = a (J-^j for 9 > 9 E j (29) 

hence, a is the surface mass density at the Einstein radius. For example, a singular 
isothermal sphere would have a = 0.5 and q= 1. Note that we do not have to specify the 
mass distribution inside the Einstein radius; it suffices to know that the mean surface 
mass density inside the Einstein radius is R = 1. Using the methods described in Chap. 
8 of Schneider et al. (1992; hereafter SEF), we find that 

/m / 2a \ / 9 \~ 2 2a / 9 \ ~ q 

and 

7(0) = k{6) - k{9) . (31) 

Thus, this family has a 'shape' parameter q and a 'strength' parameter a. 

For the second family of lens models (hereafter Family B) , we use the model described 
in Sect. 8.1.5 of SEF, which reads 

where 9 C is the core radius, and kq the central surface mass density. For 9 6 C , the 
surface mass density is nearly constant at k,q, whereas for 9 ^> 9 C it decreases like a power 
law, with a slope determined by p. The Einstein radius is given as 



9e = 9 c ^k 1 /{1 - p) -1=: 9 C VW, (33) 

and only exists if kq > 1, a well-known result for centrally-concentrated axially-symmetric 
lenses. Hence, we can eliminate 9 C in terms of #e i n (32) and again obtain a two-parameter 
model (for fixed 6>e), 

l+ P W(9/9y) 2 

K{0) = Ko [1 + W(9/9 E )^- P ' (34) 

which leaves kq and p as the two parameters. Asymptotically, for 9 — > oo, the mass 
distribution behaves like 

K^K pW p - 1 (9/9 E ) 2p - 2 . 
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Fig. 1. The radial mass profile k (upper panels), the magnification signature for the image density 
depletion, /j,~ ' 5 (middle panels), and the reduced shear g (bottom panels) as a function of angular 
separation from the cluster center, for four different combinations of the model parameters a and q, as 
indicated. Left panels are for models of Family A, right panels for Family B 



We shall choose the two free parameters in such a way that the asymptotic form of this 
model is similar to that of Family A. Hence we take a = KopW p ~ 1 and q = 2 — 2p as our 
parameters. Then, 1/W = [2a/(2 - q)] 2 / q - 1, k = (1 + W) q/2 , and therefore, 
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accordingly, 

2a - - V2 



2-q 



2a \ 2/q ( 



-1 + 



2-q - \9 



(36) 



As before the shear is calculated from (31). The condition kq > 1, or W > 0, translates 
into 

a > ^ • (37) 

For a = (2 — q)/2, the resulting lens model coincides with the corresponding one from 
Family A. 

In Fig. 1 we plot, for four different combinations of a and g, the radial mass profile 
k{9), the magnification signature for f3 = 0.5, and the reduced shear g, for both 

Families of models. The upper panels show the radial dependence of the surface mass 
density k(9) which are all fairly different over the range of angles plotted. The bottom 
panels show the radial dependence of the reduced shear; here we see, in particular for 
Family A, that models with the same amplitude but different slope have nearly the 
same curves g(9); hence, it will be quite challenging to determine the radial slope q 
for these models from the shear method. The reason for this similarity comes from an 
accidental cancellation of the slope dependencies of 1 — k and 7 in g for these models. In 
contrast to the reduced shear, the magnification signal, plotted as i$~ x = /U -0 5 in the 
middle panels, differs clearly between models with different slopes. Hence, despite the 
conclusion reached in Sect. 2 on the relative sensitivity of the shear and magnification 
methods, from the curves in Fig. 1 one might expect that the magnification effect could 
more easily distinguish between different radial slopes. 



6 Results from the likelihood analysis 

We now proceed to consider the likelihood functions for the models described in the 
previous Section, starting with the analytic results in Sect. 6.1. In Sect. 6.2 the results of 
the numerical simulations are used to substantiate our analytic treatment. We also ask 
whether it is possible to distinguish between lens Families, using each of the methods. 
In Sect. 6.3 we consider how the uncertainty in the unlensed background number density 
influences the magnification likelihood functions. 

6.1 Analytic results 

In Figs. 2 and 3, contours of constant likelihood are plotted, for (£ 7 ), (ifj), and (£tot), 
for models of Family A and B, respectively. The contours are such that one expects the 
parameters obtained by minimizing the log-likelihood from realizations of galaxy image 
observations to lie within the contours in 68.3%, 90%, 95.4%, 99%, 99.73%, and 99.99% 
percent of all cases. This expectation is based on the approximation mentioned at the 
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Fig. 2. Contours of constant average likelihood for models of Family A. The thin solid lines are contours 
of constant {£-,), the dashed contours correspond to {tfj), and the heavy solid contours to (£tot)- Contours 
are drawn for 1AI = {2.30, 4.61, 6.17, 9.21, 11.8, 18.4}, within which one expects that 68.3%, 90%, 95.4%, 
99%, 99.73%, and 99.99% respectively, of parameter estimates from realizations will be enclosed. The 
input model is described by a = 0.5, q = 1.0, and corresponds to the minimum of (£). Parameters are 
as described in the text, i.e., n M = 120arcmin~ 2 , n 7 = 30 arcmin~ 2 , (3 = 0.5, <r e = 0.2, and the inner 
and outer radii of the annulus are (a) 9 m = 0'6, # ut = 2'.0; (b) 0- m = 0'6, O ut = 4'.0; (c) 9 ln = 0'9, 
#out = 2'0; (d) 9 ln = 0'9, # ut = 4'.0. The three heavy solid curves are curves of constant (N^); the 
middle one is for (iV M ) fixed at the value obtained for the input model, the two others are for changes of 
(N„) by ±5% 

end of Sect. 4. The various panels in the two figures differ in the choice of the inner and 
outer radii of the annulus in which the data are assumed to be analyzed. 

Considering Fig. 2 first, the first thing to note is that the contours obtained from 
the shear method are very nearly parallel to the (/-axis. This is in accordance with our 
discussion of Fig. 1 at the end of the previous section where we have seen that the reduced 
shear depends only weakly on the slope q of the mass profile. Secondly, the direction 
of elongation of the likelihood contours corresponding to the magnification method is 
nicely off-set from the (/-axis; this implies that the combination of shear and magnifica- 
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Fig. 3. Same as Fig. 2, but now for models of Family B. The input model is characterized by a = 0.8, 
q = 1.1. Inner and outer radii are the same as in Fig. 2. The heavy solid line in each panel denotes the 
limit in (37); models below this line do not correspond to centrally condensed mass distributions with 

K.Q > 1 

tion information is much more powerful than each method individually in distinguishing 
between the model parameters. This fact is actually demonstrated by the likelihood 
contours corresponding to (£tot) which are much more compact than the other two sets 
of contours. Considering the upper left panel as an example, the 1-a contour of the 
magnification method closes within the boundary of the figure, whereas the 1-a contour 
of (£y) is open both to the top and the bottom; hence, from the magnification alone one 
can constrain the slope and the amplitude, whereas the shear method yields basically no 
constraint on the slope. Increasing the outer radius of the annulus, the shear method be- 
comes substantially more sensitive; in the upper right panel of Fig. 2, the areas enclosed 
by the corresponding contours for the shear and the magnification methods are about 
the same. 

In order to check whether most of the information comes from galaxy images very 
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close to the critical curves at 6 = #e = 0'.5, we have increased the inner radius of the 
annulus in the lower two panels. Although the likelihood contours widen slightly, this is 
not a dramatic effect; in particular, the likelihood contours for the combination of both 
methods remain basically unchanged. Therefore, most of the signal is not located right 
at the Einstein radius. This is good news, since there the measurement of galaxy shapes 
and fluxes may be affected by bright galaxies in the center of the cluster. 

For the models of Family B, the situation is qualitatively the same as discussed for 
Family A. Here, the advantage of increasing the outer radius of the annulus is even larger 
for the shear method than was the case for the Family A models. This result shows the 
importance of wide-field images for a quantitative analysis of cluster mass profiles (see 
Bonnet et al. 1994). As for Family A, little information is lost by increasing the inner 
radius from 1.26>e to 1.8#e- 

In order to investigate a larger part of parameter space, we consider the behaviour 
of the log-likelihood function near the minimum. Writing 

(£) (tt) » (£) (7r t ) + (Air), (An), (38) 

for parameter values close to those of the 'true' parameters, curves of constant A (£) are 
described as ellipses. The '1-a error ellipse' is described by setting 2A (£) = 2.30 = X. 
Note that from Figs. 2 and 3 one finds that the confidence contours are not always well 
approximated by ellipses; nevertheless, the consideration here should provide a qualita- 
tive handle on the expected uncertainties in parameter determinations. The matrix of 
partial derivatives in (38) has been calculated for the magnification method in (24); the 
corresponding expression for the shear method is too long to be reproduced here. The 
'1-a error ellipse' is characterized by three numbers, corresponding to the major and 
minor axes, and their orientation. In order to translate this into errors Aa and Aq for 
parameter estimates, we chose the following prescription: We construct the smallest rect- 
angle which encloses the ellipse, and take as errors the half-lengths of its sides. Denoting 
the matrix in (38) by V, one then obtains 

Aa _ / detV V? 2 I 1 

7T ~ V V&V22 + WiV V 22 detV ' 

\ , (39) 

Aq / detV V? 2 / 1 

TT ~ \Jv^Vg + T^V^iidety ' 

This graphical prescription for the error in the two parameters is certainly not optimal, 
but easy to apply; in cases where the eigendirections of V are close to 'diagonal' (as is 
the case for the magnification method in Fig. 2), this method can severely overestimate 
the error. 

The errors determined in that way are plotted in Fig. 4, for two choices of the radii 
of the annulus in which observations are assumed to be available; the small (large) outer 
radius corresponds to the left (right) panels in Fig. 4. They are plotted as a function of the 
value of a, for five different values of g, and using the shear (solid curves), magnification 
(dotted curves) and the combination of both methods (dashed curves). The most obvious 
point to note is that for the smaller outer radius of the annulus, the magnification method 
yields smaller errors in q than the shear method, for all values of a and q (see lower left 
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Fig. 4. The errors Aa and Aq, as denned in the text, plotted as a function of the input value of a. The 
lens model Family A is considered here. The left panels are for 9 m = 0'6, # ut = 2'0, and the right 
panels for ln = 0'.6, # ut = 6'0. Each panel shows three sets of five curves; the solid curves correspond 
to the shear method, i.e., are based on {£-y), the dotted curves correspond to the magnification method, 
and the dashed curves are the errors from using the combination of shear and magnification information. 
Curves within each set correspond to different input values of q, with q = 0.6,0.8, 1.0, 1.2, 1.4. For (t^) 
and (£tot), the curves vary monotonically with q, the lowest corresponding to q = 0.6, the highest to 
q = 1.4. The dependence of the errors for (£ 7 ) on q is more complicated and not necessarily monotonic; 
for small a, the uppermost curves in the left panels correspond to q = 0.6, for the upper right panel, the 
behaviour is monotonic, with Aa increasing with q. For the lower right panel, the two largest values of 
Aq correspond to the two largest values of q; the curves for the other values of q are almost coincident 



panel), whereas for the error in the amplitude a, the relative merits of both methods 
depend on a and q. For the larger outer aperture radius, the amplitude of the mass 
distribution is much better determined from the shear method, and the relative accuracy 
of the slope determination depends mainly on q: for steep profiles, q ^ 1, the shear 
method yields more accurate results, whereas the magnification method is superior for 
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flatter profiles. Except for the shear method, the accuracy increases with decreasing q 
and increasing a, or in other words, with increasing lens strength. For the smaller outer 
radius of the aperture, the flattest profiles yield the largest errors for the shear method, 
which is counter-intuitive and must be related to the accidental cancellations of effects 
in the radial profile of the reduced shear, seen in Fig. 1. 

6.2 Numerical simulations: comparison of the likelihood analysis with x 2 
statistics 

For each method, we now quantify the agreement between the distribution of 2 A (£} = 
(£) (n) — (£) (7r t ) and that expected if 2 A (£) followed a perfect \\ distribution. Ten 
thousand catalogs of lensed galaxies were generated using a Family A lens model, with 
true parameters 7r t describing the cluster, and the log-likelihood functions were minimized 
to obtain the best fitting parameters n for each realization. For each n we then calculate 
2A(£) and derive the cumulative probability distribution P(> 2A{£)), which can be 
compared with the perfect distribution P(> x|). In Fig. 5 the ratio of P(> 2A(£)) : 
P(> X2) is plotted against 2A (£) for each of the magnification, shear and combined 
methods. The deviation from a \\ distribution is very small for the magnification method, 
and slightly larger for the shear and combined methods, but still very acceptable: the 
deviation from \\ statistics is such that the recovered values of 71 for about 94.4% of the 
realizations lie within the 95.4%-confidence interval. 

The scatter of the recovered values of n in the a-q plane is also consistent with the 
ensemble-averaged log-likelihood contours, as we will demonstrate below. We can also 
ask how well each of the methods can distinguish between the Families of lens models. 
To these ends, we generated 500 catalogs of lensed galaxies using the prescription for 
numerical simulations given in Sect. 5, with a Family A lens model to describe the cluster. 
The log-likelihood functions £^, £ 7 and £tot were independently minimized under a Family 
A lens model and a Family B lens model, to obtain the best fitting parameters and their 
corresponding log-likelihood values; we refer to these as ^ M (A), £ 7 (A) and £ to t(A) where 
the minimization is performed under a Family A lens model, and ^(B), &y(B) and 4ot(B) 
where a Family B lens model is used during the minimization. By comparing the log- 
likelihood values corresponding to recovery using each of the two Families, we determined 
how frequently the minimum corresponds to recovery by Family A. 

In the upper, middle and lower panels of Fig. 6 we show our results for the magnifi- 
cation, shear and combined methods respectively. The left-hand panels show parameter 
recovery under a Family A lens model, whereas a Family B lens model was used to recover 
the best fitting parameters indicated on the right-hand panels. Contours of constant av- 
erage likelihood are marked on the left-hand panel; note that the distributions of Monte 
Carlo points are consistent with these. 

In 47% of the realizations ^(A) < ^(B) indicating that the magnification method 
cannot discriminate between Family A and Family B. The middle panel of Fig. 6 shows 
that the shear method performs slightly better: in 64% of the realizations £ 1 (A)< £ 7 (B). 
The results for recovery using the combined method are shown in the bottom panel; 
^tot(A)< 4ot(B) in 63% of the realizations. 

We repeated the same procedure, generating catalogs with a Family B lens model, 
and performing the minimization of i^, £ 1 and £ to t under a Family B lens model and 
under a Family A lens model; our results are shown in Fig. 7. Contours of constant 
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Fig. 5. For each of the magnification (lower curve), shear (upper curve) and combined methods (dashed 
middle curve), we show the ratio of P(> 2A (£)) to that expected if 2A (£) followed a x 2 distribution 
versus 2 A (£}. On the top edge of the plot the 68.3%-, 90%- and 95.4%-confidence intervals are marked 

average likelihood are marked on the left-hand panels and the Monte Carlo points are 
consistent with these. 

Again, the shear method (middle panel) fares a little better than the magnification 
method (top panel) in discriminating between the Families: in 57% of the realizations, 
£ )Jl (B)< ^ M (A) and in 73% of the realizations £ 7 (B)< £ 7 (A). The combined method (bot- 
tom panel) gives 4ot(B)< 4ot(A) in 78% of the realizations. 

6.3 Uncertainties in the unlensed number counts 

Up to now we have assumed that the number density of background galaxies at the flux 
threshold of the observations is known. In fact, the shape of the likelihood contours 
in Figs. 2 and 3 are such that they closely follow the curve of constant expected total 
galaxy number in the annulus, shown as the middle of the heavy dashed curves in those 
figures which are curves of constant (A M ) = (N fl ) t . This implies that to first order, the 
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Fig. 6. A Family A input model (a = 0.8, q = 1.3) was used to generate lensed galaxy catalogs. For each 
of the magnification (top panel), shear (middle panel) and combined methods (bottom panel) the left- 
hand (right-hand) panel shows parameter recovery under a Family A (Family B) model and the crosses 
(triangles) indicate where the log-likelihood function is lowest when the minimization is performed using 
a Family A (Family B) lens model. On the left-hand panels we mark contours of constant A {£) with 
levels the same as in Fig. 2 
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magnification information is provided by this total number. As indicated by the other 
two heavy solid curves in these figures, a change of the number by ±5% changes the curve 
of constant (N^) quite drastically. Given that a high precision of the number density 
of very faint galaxies may be difficult to achieve, we shall now consider two different 
assumptions: first, we shall assume that the unlensed number density is not determined 
from any calibration frame, but is obtained from the present data set itself. After that, 
we shall consider the (more realistic) assumption that the unlensed number counts is 
known, but with some uncertainty. To this end, we reinsert the density-dependent term, 
which we discarded in (11), into the log-likelihood function, which then reads 

e^ = n^i + (i-p)J2 ln M0*) - % ln ?v , ( 4 o) 

where we have defined 

r d 2 eHe)f- 1 . (4i) 



This expression can now be minimized with respect to the unknown density n M , which 
yields, not unexpectedly, 

n n = — • 

Inserting this value into (40) yields the new log-likelihood function, up to a model- 
independent constant, 

4 = (1 - (3) X>M*i) + N„ In / . (42) 
Following the techniques applied in Sect. 4, the ensemble average of £ M can be calculated, 



M = n 



d 2 9 M*)]"- 1 ) In/ 
+ rv(l-/3) J d 2 6 [MO)] ' 1 In 



(43) 



where now n M is the true number density of galaxies. By differentiation with respect to 
7Tj, one can easily show that (^ixj has a minimum when n = n t . 

In Fig. 8 we have plotted contours of constant likelihood of ^m}' f° r two vames °f 

the outer radius of the annulus, together with the corresponding contours of (£ 7 ), 

and (itotj '■= (i-y F rom the figure it becomes immediately clear that dropping 

the assumption about the knowledge of the unlensed number density leads to a drastic 
loss of information from the magnification method. The thin dashed contours in Fig. 8 

which correspond to (&i^J are much wider than the solid contours corresponding to 

In addition, whereas the contours of (£^) were seen to be substantially misaligned 
relative to those of (£ 7 ), so that the combination of shear and magnification information 
substantially decreases the error region in parameter space relative to any single one 

of these methods, the contours of ( j are basically parallel to those of (£ 7 ) (shown 
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Fig. 7. A Family B input model (a = 0.8, q = 1.1) was used to generate lensed galaxy catalogs. For each 
of the magnification (top panel), shear (middle panel) and combined methods (bottom panel) the left- 
hand (right-hand) panel shows parameter recovery under a Family B (Family A) model and the crosses 
(triangles) indicate where the log-likelihood function is lowest when the minimization is performed using 
a Family B (Family A) model. On the left-hand panels we mark contours of constant A {£} with levels 
the same as in Fig. 2 
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Fig. 8. Similarly to Figs. 2 and 3, contours of constant A {£} are plotted; the contour levels are the same 
as described in Fig. 2. The inner radius of the annulus is 6 m = 0'6, the outer radius is 6 ou t = 4' (upper 
panels) and 15' (lower panels); left (right) panels correspond to models of Family A (Family B). The 

thin dashed contour correspond to (^/n), the thin solid contours to (l/j,), the heavy dashed contours to 

(£-y), and the heavy solid contours to ^tot^) := +£fi)- The input model is a = 0.5, q = 1.0 for 
Family A, and a = 0.8, q = 1.1 for Family B 

as heavy dashed curves in Fig. 8). Therefore, the combination of both methods, in the 
absence of knowledge of the unlensed number density, yields only slightly more accurate 
parameter estimates than the shear method alone. 

To understand this point, one needs to realize why the magnification information 
turned out to be so powerful when is known. As we have seen in Fig. 1, the reduced 
shear of models with different q is very similar, owing to the mass-sheet degeneracy 
discussed in the introduction. Knowing n M , the mass-sheet degeneracy is broken, which 
leads to the sensitivity of the magnification to the slope of the density profile. However, 
if is unknown, the mass-sheet degeneracy remains, and the discussion of the relative 
merits of shear and magnification methods in Sect. 2 prevails. 
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We now make the more realistic assumption that the unlensed number density is 
known approximately. Suppose one knows the value of with a fractional accuracy of 
77; one can then supplement the likelihood fucntion with a prior, taken to be a Gaussian 
in n M , with mean n M , and dispersion a n = rjn^. Hence, the corresponding log-likelihood 
function becomes 

= n„ I + (1 - ff) ^ In Mi) ~ In n, + ( ^^f ) ■ (44) 
As before, this can be minimized with respct to n M , yielding 

^ = 5 (1 - + ^(l-^) 2 +^iV At . (45) 

Inserting this value for n M into (44) yields the new log-likelihood function, which we shall 

denote as £^ . Taking the ensemble average in the present case is slightly more difficult, 
owing to the occurrence of iV M in the square-root. However, if we neglect the Poisson- 
fluctuation in the ensemble averaging in (17), which is a very good approximation when 
{Nfj) t is large, the ensemble average can be obtained as 

(W) = M I + n„(l -f3) j d 2 9 rf" 1 In/, - (N,) t In (n M > + . (46) 

where 

^ = \ (1 - ^M) + yj^-rf^if + n* Wt ■ (47) 

Although this expression looks fairly complicated, it is easy to show that (ift^ has 

a minimum when 7r takes the value of the true model. For 77 — > 0, this log-likelihood 
reduces, up to additive constants, to the one for which the unlensed number density was 
assumed to be known, i.e., eq.(ll), whereas for large 77, it approaches (43). 

In Fig. 9 we have plotted the 90%-confidence contours as obtained from (^^^, 

using models of Family A, for various outer radii of the data annulus. Four values of 
the fractional uncertainty 77 in the number density have been chosen, ranging from 77 = 
- for which the contours agree with the corresponding ones in Fig. 2 - to 77 = 0.06. In 
addition, the 90%-confidence contour for the shear and those for the combined method 
are plotted. As expected, by increasing 77, the confidence region increases, and the relative 
increase is larger, the larger 9 out - The reason for this is that for the small annulus, the 
relative change of the number density by magnification is much larger than the fractional 
uncertainty, whereas this is no longer true at larger radii. In fact, as can be seen from the 
lower right panel in Fig. 9, the confidence contours seem to approach an asymptotic form 
quickly for large annuli, implying that the prior information becomes irrelevant unless 77 
is very small, since the unlensed number density is then more accurately determined from 
the data set itself. Thus, the outer-most contour in the lower-right panel in Fig. 9 is very 
similar to the corresponding one in Fig. 8, where no prior information about the number 
density was assumed. Thus, the larger the region from which data are used to determine 
the mass profile parameters, the more accurate the unlensed number density needs to 
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Fig. 9. For the case that the unlensed number counts are known with a fractional accuracy 77, the light 
dashed contours plot the 90%-confidence contours of parameter estimates as obtained from the ensemble- 
averaged log-likelihood function (46). The model a = 1/2, q = 1 of Family A has been taken as input. 
The inner radius of the annulus in which data are assumed to be is #j n = 0'.6, and the outer radius is 
different in each panel. The solid contour is the 90%-confidence contour as obtained from the shear, and 
the heavy dashed curves from combining shear and magnification constraints. The sets of four dashed 
curves are for rj = 0, 0.02, 0.04, and 0.06. 



be known, for a given accuracy relative to that obtainable from the shear method. In 
particular, we conclude that the magnification method considered here seems to be most 
powerful for data sets which do not extend far from the cluster center, whereas for wide- 
field imaging data, the shear method is superior. 

Besides observational difficulties to determine the unlensed number density with a 
given selection function from calibration data, the anglar correlation function of galaxies 
provides a fundamental limit down to which the unlensed number density can be de- 
termined within any given region. To this end, we have calculated the fluctuations of 
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the number of galaxies in a circular aperture of radius 9, using a two-point correlation 
function of the form u;(i?) = A# -0 8 . As a result, we find that 5N/N oc 9~ 0A , a slowly 
decreasing function of the radius. Unfortunately, at the faint magnitude limits where 
the magnification method is applied, no very accurate determination of the angular two- 
point correlation function is available; the best measurements has been carried out in 
the Hubble Deep Field (e.g., Villumsen et al. 1997, and references therein). Whereas 
at B ~ 28 a clear signal in the correlation function is seen on scales below ~ 1 arcsec, 
the signal is too weak, and the solid angle of the HDF too small to determine whether 
the power-law behaviour of u extends beyond arcsecond scale. If one extrapolates the 
measurements of Villumsen et al. (1997) to larger angular scales, then 5N/N ~ 0.1 on 
an angular scale of ~ 4'. In that case, the results in Fig. 9 suggest that the value of the 
magnification method is rather limited. More likely, though, the preceeding value is a 
vast overestimate of the number count fluctuations due to angular clustering at these 
magnitudes. 

7 Discussion and conclusions 

Using maximum likelihood techniques, we have investigated the accuracy with which 
weak gravitational lensing can recover the radial density profile of clusters of galaxies. 
To keep the treatment as simple as possible, we have assumed that the Einstein radius 
of the cluster is known from strong lensing features, and that only the weak lensing 
information from outside the Einstein radius is used; we also confined the consideration 
to axially-symmetric mass distributions only. We have compared the shear method with 
the magnification method in form of the number depletion. Although from a simple 
analytical consideration the latter method is expected to yield less accurate results, the 
mass sheet degeneracy causes the reduced shear to be nearly degenerate along a one- 
dimensional subset of the two-parameter mass models considered here. In some cases 
it then turns out that the magnification method can constrain the parameters of the 
mass model more accurately than the shear method. However, this is true only as 
long as the unlensed number density of galaxy images is assumed to be known quite 
accurately. In the absence of this knowledge, the magnification method is inferior to the 
shear method. Since the number density of galaxies obtained from 'blank field' imaging 
depends not only on the flux threshold (which requires accurate photometry) , but also on 
the seeing and sky brightness, this accurate number density calibration may be difficult 
to achieve. The advantage of the shear method is that no external calibration is required 
(the external calibration is replaced by the assumption of random intrinsic orientations). 
In addition, the angular correlation of faint galaxies may present a fundamental obstacle 
for an accurate determination of the unlensed counts in a given area, but too little is 
known about the angular correlation function at flux limits at which the magnification 
method is applied to make quantitative estimates of this effect at present. 

Without accurate external calibration, our results show that large areas around 
the cluster need to be mapped, since the external calibration is effectively replaced by 
calibration from the outer region of the data field. In order for this calibration to be 
accurate the image quality needs to be homogeneous over the wide field. Of course, since 
the measurement of shear requires superb observing conditions, in particular with regards 
to the seeing, the number count magnification method may be applied to photometric 
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data set which cannot be used for the shear method. 

Although we have not investigated the magnification effect based on the change of 
image size at constant surface brightness (Bartelmann & Narayan 1995), the previous 
remarks on external calibration apply there as well. 

External calibration may be possible using HST observations where the observing 
conditions are much more stable than from the ground. Hence, the claim by Broadhurst 
(1999) that the mass profile of the cluster C10024+17 can be better constrained with the 
magnification method than with the shear method is not contradicted by our results here; 
in fact, if the external calibration of the number density is possible for these HST images, 
the results shown in Figs. 2 and 3 indicates a higher accuracy of the shape parameter 
than possible with the shear method. 

We would like to thank Matthias Bartelmann and Tom Broadhurst for several fruitful 
discussions, and also the former for carefully reading the manuscript. This work was 
supported by the TMR Network "Gravitational Lensing: New Constraints on Cosmology 
and the Distribution of Dark Matter" of the EC under contract No. ERBFMRX-CT97- 
0172 and the "Sonderforschungsbereich 375-95 fiir Astro-Teilchenphysik" der Deutschen 
For schungsgemeinschaft . 



Appendix 

We derive here the dispersion a of the image ellipticity distribution. As in (2) let p e (e s ) 
be the probability density of the intrinsic source ellipticities. Without loss of generality, 
we choose the reduced shear g to be real, and write e s = xe lip . Then, the transformation 
(1) reads 



2 gx + cosp 

e 1 =g + (l-g )x — — ■ — ^ 

1 + 2gx cos <f + g z x z 

€ 2 = (l-g 2 )x 



X 6 p e (x) 2 

x A g z 



(Al) 



1 + 2gx cos ip + g 2 x 2 
Then, the dispersion of the observed ellipticities in the two directions becomes 

p 1 /»27T 

cr| = / dxxp(x) / d<p (e 1 -g) 2 
Jo Jo 

=a-^jf *> ( 1+2 ::r;:^ ) 2 ^ 

= (l- 9 2 )V^ 1 d I ^ 
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